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ABSTRACT 

Nebular-phase spectra of SN2006aj, which was discovered in coincidence with 
X-ray flash 060218, were obtained with Keck in 2006 July and the Very Large 
Telescope in 2006 September. At the latter epoch spectropolarimetry was also 



attempted, yielding an upper limit of ~2% for the polarization. The spectra show 
strong emission lines of [O i] and Mg i] , as expected from a Type Ic supernova, but 



weak Call lines. The [Fell] lines that were strong in the spectra of SN1998bw 
are much weaker in SN2006aj, consistent with the lower luminosity of this SN. 
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Based on observations made with ESO Telescopes at the Paranal and La Silla Observatories under 
O program 277.D-5039(A). 
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The outer velocity of the line-emitting ejecta is ~ 8000 km s -1 in July and ~ 
7400kms _1 in September, consistent with the relatively low kinetic energy of 
expansion of SN2006aj. All emission lines have similar width, and the profiles 
are symmetric, indicating that no major asymmetries are present in the ejecta 
at the velocities sampled by the nebular lines (v < 8000 km s -1 ), except perhaps 
in the innermost part. The spectra were modelled with a non-LTE code. The 
mass of 56 Ni required to power the emission spectrum is ~ 0.20 M Q , in excellent 
agreement with the results of early light curve modelling. The oxygen mass is 
~ 1.5 Mq, again much less than in SN 1998bw but larger by ~ 0.7 M & than the 
value derived from the early-time modelling. The total ejected mass is ~ 2 M & 
below 8000 km s -1 . This confirms that SN2006aj was only slightly more massive 
and energetic than the prototypical Type Ic SN 19941, but also indicates the 
presence of a dense inner core, containing ~ 1 M & of mostly oxygen and carbon. 
The presence of such a core is inferred for all broad-lined SNe Ic. This core may 
have the form of an equatorial oxygen-dominated region, but it is too deep to 
affect the early light curve and too small to affect the late polarization spectrum. 

Subject headings: supernovae: general — supernovae: individual (SN2006aj) - 
nucleosynthesis — gamma rays: bursts 



1. Introduction 



The connection between some gamma-ray bursts (GRBs) and supernovae (SNe) is now 



well established (see IWoosley fc Bloom! 120061 for a recent review). The first observational 



clue for a connecti on came with the discovery of a brigh t SN coincident in time and space 
with GRB 980425 JGalama et alill998l : Irian et al.lboooh The SN was eventually classified 
as a Type Ic (for a definition of the class, see iFilippenkd 119971 ). indicating an origin from a 
massive star that lost both the H and He envelopes prior to core collapse. SN 1998bw was an 
exceptional SN Ic in being very luminous and showing a broad-lined spectrum , indicative of a 



large kinetic energy of expansion. Such objec ts have been called hype rnovae (jlwamoto et al. 



19981 ) or broad-lined SNe Ic (SNe Ic-BL; see lWoosley fc Bloom! 120061 ) . Modelling confirmed 
the large kinetic ener gy, ^ (10 — 30) x 10 51 erg, and showed that the ejected mass was also 
quite large (~ 10 M g . llwamoto et al.lll9981 ). pointing to a progenitor star with zero-age main 



sequence (ZAMS) mass ~ 40 M . 

Later, other nearby long-duration soft-spectrum GRBs, in the same high-energy class as 
GRB 980425, were also seen to accompany hypernovae/SNe Ic-BL, confirming the connection 
( jStanek et al.ll2003l ; iMalesani et al.ll2004l ). The SNe had properties very similar to those of 
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SN 199 8bw, although the properties of the GRBs were quite different (IMazzali et al.l 12003 
2006al ). In a popular model, both the SN and the GRB a re produced when the rapidl y 
rotating core of a massive star collapses to form a black hole (IMacFadyen fe Woosleylll999l ). 
The coincidence of a very aspherical event, the GRB, and a SN motivates the search for 
signatures of asphericity in the SN as well. 

In the case of SN 1998bw, evidence of asphericity came from the nebular spectra. These 
showed strong emission lines of both [Ol] AA6300, 6363 and [Fell] (a blend near 5200 A), as 
could be expected from a bright SN Ic, but with the peculiarity that even after accounting 
for blending the [Fe n] lines were broader than the [O i] line, indicating a higher expansion 
velocity of Fe than of O. This was explained in the context of an axisymmetric explosion, 
where the Fe produced by the decay of 56 Ni is ejected at high velocities near the direction 
of the GRB jet while oxygen, which is mostly not a product of the explosi on but rather 



left o ve r from the pr o genito r, is ejected at lower velocities nearer the equator (IMazzali et al. 



20011 ). iMaeda et al.l (120021 ) supported this scenario with two-dimensional explosion models 
and three-dimensional synthetic nebul ar spectra, and sug gested that our viewing angle was 
~ 15 — 30° from the jet axis fse e also IMaeda et al.l 120061). This may partly be responsible 
for the weakness of GRB980425 ([Ramirez-Ruiz et al.ll2005l ). 



In the case of SN 2003dh/GRB030329, polarization levels of ~ 2% were observed in 
the fi rst 2-3 days after the GRB, but they were probably measurements of the GRB after- 
glow (IKawabata et al.ll2003l ; iGreiner et al.ll2003l . and references therein), dominated by the 
synchrotron process which imposes partial polarization. Later polarization measurements, 
which should not have been affected by the afterglow, suggested that the polarization was 
small (P < 1%), al though the alignment of the polarization angle suggests that the polar- 
ization is non-zero (IKawabata et al.ll2003l ). This could be consistent with a jet-like explosion 
viewed very close to the jet the strength of GRB 030329 may suggest. 

A later confirmation that the ejecta of bright SNe Ic can be very aspherical came from 
the detection of the d ouble-peaked emission profile of the [O i] AA6300, 6363 line in SN 2003jd 
(IMazzali et al.ll2005l ). The profile could be reproduced assuming that we observed an explo- 
sion similar to that of SN 1998bw from close to the equatorial plane. A GRB may or may not 
(jSoderberg et al.l l2006ah have been produced, and if it was it may not have been detected 
because it was not pointing toward the Earth. 



Asphericity may indeed be a common feature of core-collapse SNe ([Leonard et al.ll2006l ). 
and it is probably stronger i n the carbon-oxygen cores, since envelope-stripped SNe show a 
highe r degree of polarization (jWang et al. l ll996l : lFilippenko fe Leonardlbooi [Leonard fe Filippenkol 
2005) . Spectropolar i metry of the non-G RB, broad-lined SN Ic 2002ap ( Kawabata et al. 20021 ; 
Leonard et al.ll2002l ; IWang et al.ll2002l ) revealed polarization of up to ~ 2% in the O I A7774 
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line and in the Call IR triplet at early times, suggesting a deviation from sphericity of 
~ 10 - 20%. 



X-ray flashes (XRFs) are a lower-energy and softer subclass of long-duration GRBs 



Heise et al 



( jFynbo et al. 



2001). Previous work had indicated a connection between XRFs and SNe 



2004 ; iTominaga et al.l 12004 ; ISoderberg et al.l 120051 ) but high-quality spectra 



of the SNe had not been obtained until the Type Ic SN200 6aj was discovered coincident 
with XRF 060218, at a redshift z = 0.03342 Jpian et al.ll2006h . Although SN2006aj showed 
a broad-lined spectrum, this was not as broad as those of SNe 1998bw and 2003dh. There- 
fore, SN2006aj could not, prima facie, be considered an extremely energet ic SN from early 



optica l spectroscopy and photometry. Radio observations confirmed this (ISoderberg et al. 
2006bh . 



SN2006aj was modelled (IMazzali et al.l l2006bl ) to have synthesized only ~ 0.2 M & of 
56 Ni, and to have ejecta only slightly more massive (~ 2 M ) and energetic (E pa 2 x 10 51 erg) 
than "normal" SNe Ic0 These parameters a r e more appropriate for the collapse of a star 
of ~ 20 M to a neutron star. IMazzali et al.l (J2006bJ) suggested that magnetic activity on 
the surface of the nascent neutron star increased the energy of the explosion and caused the 
XRF. In this case, the explosion may not be highly as p herical . The absence of a jet break 
and the behavior of the radio flux led ISoderberg et al.l (l2006bl ) to estimate a "jet" opening 
angle of > 50°. 



Large polarization (~ 4% at 3-5 days) was observed in SN2006aj (IGorosabel et al. 



20061 ). As in the case of SN2003dh, this may have been due to the relativistic jet and the 
XRF afterg low synchrotron contin uum, but as the data have a later epoch than those of 
SN2003dh (IKawabata et al.l 120031 ). and the afterglow was not as strong as in the case of 
SN2003dh, an aspherical fast outflow from the o utermost SN ejecta may also have been 
detected. At later epochs, IGorosabel et al.l (120061 ) find that the polarization level dropped 
to ~ 1.4%, which they attribute to the host galaxy. 

Obviously, it is important to verify these suggestions using nebular spectroscopy. We 
observed SN 2006aj with the Keck I telescope and the Very Large Telescope ( VLT) when it 
became visible again after solar occultation. In the following, we discuss the observations 
(§2), present the observational results (§3) and the models of the nebular spectra (§4), and 
discuss our findings (§5). 



14i 



i We use the term "normal" to refer to SNe like SN 19941 ( Sauer et al. 20061 and references therein), i.e. 
SNe with explosion energy « 10 51 erg, although it is by no means clear that this can be called a normal 
event. 
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2. Observations and Spectral Analysis 



On 2006 July 26.61 (UT dates are used throughout this paper), SN2006aj was observed 
with the 10-m Keck I telesco pe on Mauna Ke a, Hawaii, equipped with the Low Resolution 
Imaging Spectrograph (LRIS; lOke et al.lll995l ) with a D560 dichroic to separate the blue and 
red light paths. The seeing was about 0.8", and a long slit of width 1" was aligned along a 
position angle of —82°. Although the parallactic angle (Filippenko 1982) at that time was 
—52°, the atmospheric dispersion was not very significant, given the relatively low airmass 
(1.3). The blue-side data were taken with the 600/4000 grism, giving a spectral resolution 
of 3.5 A, while the red-side data were taken with the 400/8500 grating, giving a spectral 
resolution of 6.4 A. Two 600 s exposures were obtained. Standard CCD processi ng was 
accom plished with IRAFr^l. The data were extracted using the optimal algorithm of Home 
(119861 ) . Wavelength calibration was obtained from internal arc lamps and slight adjustments 
were derived from night-sky lines. The spectrophotometric standard stars BD+284211 and 
BD+ 174708 (observed at the parallactic angle), as w ell as our own IDL ro u tines, were used 



for fl ux calibration and telluric absorption removal fjWade Hornd Il988l ; iMatheson et al. 

2oooh . 



On 2006 September 19.33 and 20.36, SN2006aj was observed with the 8-m VLT Kueyen 
of the European Southern Obse rvatory (ESQ) at Cerro Paranal, equipped with the FOcal 
Reducer Spectrograph (FORS1. lAppenzellerlll998l ) in polarimetric mode. The observations 
were performed with the 300V grism and a 1.1" wide slit, aligned along the parallactic angle, 
giving a dispersion of ~2.6A pixel -1 and a resolution of ~llA at 5500 A (full width at 
half-maximum intensity). The wavelength coverage achieved with this setup is 3300-8600 A. 
At both epochs the SN was observed at four half-wave plate position angles (0°, 22.5°, 45°, 
and 67.5°), for a total integration time of 7200 s. The seeing was 1.2" on September 19 and 
1.5" on September 20. 

The data were bias-subtracted, flattened, and wavelength-calibrated using standard 
tasks within IRAF. Stokes parameters, linear polarization degree, and position angle were 
computed using specific routines written by the authors. Finally, polarizati on bias correction 
and er ror estimates were performed following the prescriptions described by lPatat &: Romaniello 
rt2006[ ) . while the half- wave plate zeropoint angle chromatism was corrected using FORS1 
tabulated datfl In order to reduced noise, the two data sets have been combined and the 
final Stokes parameters rebinned to 75 A bins (29 pixels). The signal-to-noise ratio (S/N) 



15 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Asso- 
ciation of Universities for Research in Astronomy, under contract with the National Science Foundation. 

16 See |http : //www . eso . org/ instr ument s/f or s/ inst/pola.html 
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in the binned flux spectrum is ~100 in the [Ol] AA6300, 6363 region, corresponding to an 
expected root-mean square uncertainty in the polarization degree of ~ 0.7%. 

Photometric calibration of the spectra was achieved by observing spectrophotometric 
standard stars with full polarimetric optics inserted. Instrumental polarization and position 
angle offset were checked by observing polarized and unpolarized standard stars, within the 
FORS1 calibration plan. 

The total- flux spectrum was obtained combining all exposures, giving S/N « 30 in the 
continuum at 6000 A. At this advanced nebular phase, the only relevant SN signal is expected 
to be associated with the nebular emission lines. Therefore, we attributed all observed 
continuum to the host galaxy. The continuum was evaluated in four spectral regions (near 
4400 A, i.e., immediately blueward of the Mgl] line; at 6200 A and 7020 A, i.e., on either 
side of the [Ol] line; and at 8170 A), and was then interpolated at all wavelengths of our 
spectral range with a spline and subtracted, along with the host-galaxy emission lines, from 
the observed spectra. We thus obtained the pure SN emission line spectrum. 

In order to compute line intensities the spectral flux was converted to physical units. 
The SN spectrum is superposed on the host- galaxy continuum Galaxy photometry, however. 



includes narrow host-galaxy emission lines (jCool et al.ll2006l ; iModjaz et al.ll2006l ; iPian et al. 



20061 ; ISollerman et al.ll2006l ). The flux level of the SN spectra is however best established by 
comparing to the galaxy continuum. Therefore, first we estimated from the observed spectra 
the fractional contribution of the host-galaxy emission lines to the galaxy photometry in 
each filter, and hence determined the photometric level of the galaxy spectral continuum. 
This allowed us to evaluate the flux level of the observed and subtracted SN spectra, and 
ultimately to estimate the intensities of the SN nebular emission lines. As the host galaxy 
of SN2006aj is very small and compact, the fraction of galaxy and SN light possibly not 
included in the slit during spectroscopy should be similar, ensuring that our method does 
not substantially underestimate the galaxy contribution. Figure 1 shows the observed VLT 
spectrum, the fiducial galaxy spectrum, the observed spectrum recalibrated according to the 
galaxy photometry, as explained above, and the subtracted SN spectrum, calibrated and 
cleaned of cosmic rays. 



3. Results 

Close inspection of the VLT-FORS1 polarization data shows no detection of linear 
polarization, either in the continuum or in the most prominent emission line, [Ol] AA6300, 
6363; only a conservative 3a upper limit of 2% across the whole wavelength range could be 
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set. Considering that a polarization level of P ~ 1.4% may be attributed to the host galaxy 
jGorosabel et al-lbood ). and that interstellar polarization may contribute a similar amount 
[Pis,max — 0.9 x E(B — V)], this is consistent with unpolarized SN flux, and with a global 
ellipt icity of the inner SN ejecta (i.e. at v < 8000kms _1 ) of no more than 10% (iHoflich 
19911 ). This is not expected to be observable in the profiles of the nebular emission lines. 



We derived SN magnitudes B « 22.8, V « 21.6, R « 21.0, and I « 20.9 for the Keck 
data (rest-frame epoch 153 days), and B « 23.4, V « 22.6, and « 21.8 for the VLT data 
(rest-frame epoch 206 days). We assign an uncertainty of ~ 0.1 mag in the R, V, and I 
bands, and 0.2 mag in the B band. The average light curve decline rate in this period is then 
~ 1.7 ±0.7, 1.5 ±0.4, and 1.2 ± 0.3 mag (100 days) -1 in the V, R, and I bands, respectively. 
This is consistent with expectation for a late-time light curve powered by gamma rays from 
56 Co decay in r elatively transpar ent SN ejecta and is also consistent with other SNe Ibc such 
as SN1998bw JPatat et alJboOlh and SN2002ap jFolev et alil2003l : iTomita et alJbooj ). 



The Keck spectrum, obtained on 2006 July 26, has a rest-frame epoch of 153 days. 
Figure 2 shows a comparison with the SN 1998bw spectrum obtained on 1998 September 12, 
at an epoch of 139 days after the explo sion, which is ass umed to have occurred on 25 April 
1998, in coincidence with GRB 980425 (IPian et all 120001 ). The two spectra are very similar, 
but obviously the [Ol] line is stronger relative to all other lines in SN2006aj. 

The VLT spectrum, obtained on 2006 Sept 19-20, has a rest-frame epoch of 206 days. 
It is shown in Figure 3 compared to the spectrum of SN 1998bw obtained on 1998 November 
26, at a rest-frame epoch of 216 days after the explosion. Again, [Ol] AA6300, 6363 is by 
far the strongest emission line, and is stronger in SN 2006aj than in SN 1998bw relative to 
the other lines, indicating that oxygen dominates the composition of the ejecta. The next 
strongest lines are near 4600 A, corresponding to Mgi] A4571 A, and a line near 7400 A. This 



is too red to be Call] AA7291,7324, and [Nin] A7380 is the best candidate flMaeda et al. 



20071 ). The wavelength of the emission near 8700 A in the Keck spectrum suggests that the 
Call IR triplet is not the only contributor, and that [Cl] A8727 is also strong. The [Fell] 
lines are relatively weak, but still form an emission feature near 5300 A. The strongest lines 
in that complex have rest wavelengths 5159, 5262, 5273, and 5334 A. [Fein] lines are very 
weak, as in all SN Ic spectra, indicating low temperature and significant clumping. 

The profile of the emission lines in the nebular spectrum of SN2006aj is not very sharp, 
indicating that the nebula is to a good approximation spherically symmetric. Figure 4 shows 
a comparison of the [Ol] AA6300, 6363 line of SNe 1998bw, 2006aj, and 2003jd. The line in 
SN2006aj is broader than the other two by ~ 2000 km s , although the ejecta velicity of 
SN2006aj was significantly lower than that of SN 1998bw. The line of SN 1998bw is quite 
sharp, which was interpreted as a disk-like distribution of slow-moving oxygen viewed from 
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a near-polar direction (IMaeda et al.ll2002l ). In SN2003jd the line has a width similar to that 
of SN 1998bw, but t he double-peak profile suggests that we are viewing the oxygen-rich disk 
from near its plane (IMazzali et al.ll2006bl ). 



Models 



We modelled the nebular sp ectra using a no n-LT E code (jMazzali et al.l 120011) ba sed on 
the approximations discussed in lAxelrodl (119801 ) and iRuiz-Lapuente &: Lucyl (ll992l ). The 
deposition of the gamma rays emitted in the decay of 56 Co to 56 Fe is computed using a gray 
opacity, while the positrons that are also produced are assumed to deposit their energy in 
situ. Collisional heating by the fast particles produced by the deposition of the gamma rays 
and the positrons is computed, and is balanced by cooling via line emission. The emission 
lines depend on the composition; in the case of a SN Ic, oxygen dominates the cooling. 

For both spectra the spherical approximation seems reasonable, and so we use the code 
in its simplest version, assuming that the emitting volume is spherical and homogeneous. 
This allows us to test the basic geometric properties of the ejecta of SN2006aj by verifying 
any deviation of the nebular lines from a paraboli c profile. We ado pted a distance of 140 Mpc 
and a reddening correction of E B _ V = 0.13 mag (IPian et al.ll2006l ). which takes into account 
both Galactic and intrinsic absorption. 

The main values of the fits are shown in Table 1. The outer velocity of the model for the 
Keck spectrum is 8000kms _1 , while for the later VLT spectrum this value is 7400kms _1 . 
A slight decrease of the outer velocity may be expected as the expansion makes the density 
decrease with time, such that outer regions may progressively become too thin for collisional 
excitation processes to work efficiently. 

The Keck spectrum (Figure 5) requires a 56 Ni mass of 0.20 ± 0.01 M , in excellent 



agreement with the early light-curve model (IMazzali et al.l l2006bl ). and an oxygen mass 
of 1.50 ± 0.15 M . Uncertainties are estimated from parameter combinations that yield 
acceptable fits. The line near 7400 A is mostly due to [Nill] A7380. Given the late epoch 
of the spectra, stable 58 Ni must be responsible for the emission, which can be reproduced 
assuming a 58 Ni mass of 0.02 M Q . The consequ ences of this on the p roperties of the progenitor 



star are discussed in an accompanying paper (IMaeda et al.l 120071 ). Other lines contributing 
to the emission are [Fell] AA7155, 7453 and [Coll] A7541. The Call] AA7291, 7324 lines 
are very weak. Similarly, the emission near 8600 A is mostly due to [Cl] A8727, while the 
Ca II IR triplet is very weak. The line can be reproduced only if we assume the presence of 
~ 0.3 M & of carbon. This is an interesting result, as in the early phase no carbon was visible 
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in the spectrum (IMazzali et al.ll2006bl ). The low calcium abundance is also interesting, as 
it may be a clue of the nucleosynthesis in the progenitor and in the explosion. The Ca/O 
fraction is ~ 10~ 4 by mass, whereas the corresponding value for SN 1998bw was ~ 10~ 2 . The 
total mass enclosed within the outer velocity of 8000 km s -1 is 2.07 ± 0.20 M Q . 

The 56 Ni mass estimated for the VLT spectrum is 0.19±0.01 M G and the oxygen mass is 
1.42 ± 0.15 Mq (Figure 6). The total mass enclosed within the outer velocity of 7400 km s -1 
is 1.94 ± 0.25 Mq. The slight reduction in ejected mass derived from the later VLT spec- 
trum with respect to the earlier Keck spectrum is a consequence of the slightly narrower 
emission lines, and is consistent with the model used to fit the early light curve and spectra 
( IMazzali et al.ll2006bl ). However, while the mass of 56 Ni derived from the nebular models 
(« O.2OM ) is in excellent agreement with the value estimated from the early light-curve 
evolution, the ejected mass, and in particular the masses of oxygen and carbon, are larger. 



The last early-time spectrum, obtained on 2006 March 10 (IPian et al.ll2006l ). was mod- 
elled for a photospheric velocity of 10,000 km s _1 , so there is no overlap between the region 
sampled by the early data (above 10,000 km s -1 ) and that sampled by the nebular spec- 
tra (below 8000 km s -1 ), except that the early light curve depends on the density of the 
innermost reg i ons. T he ejected mass predicted below 8000 km s -1 by the model used by 
Mazzali et al.l (l2006bl ) is ~ 1M Q , so the nebular results suggest the presence of a dense 
inner region containing an additional ~ 1M mostly composed of oxygen (~ 0.7 Mq) and 
carbon (~ 0.3 M Q ), signs of which were not visible in the early-time spectra. The presence 
of an inner high-density core is typical of other SNe Ic (e.g., SN 1998bw, SN 2002ap), as 
discussed in the next section. 

The only deviation from a parabolic profile may be seen in the core of the [O i] line. 
This may be due to the assumption of a constant density, or to the presence of an oxygen- 
dominated high-density inner core. The profile suggests that the oxygen abundance is drop- 
ping at velocities between 2000 and 4000 km s -1 , but is higher again at the lowest velocities 
(v < 2000 kms" 1 ). Between 2000 and 4000kms _1 intermediate-mass elements may dom- 
inate, but the high flux at the lowest velocities may indicate the presence of a disk-like 
oxygen-rich region. Further data and more detailed modelling will clarify this. 



5. Discussion 

The nebular spectra of SN2006aj are dominated by strong [Ol] AA6300, 6363 emission, 
as expected for a SN Ic. The masses derived from the modelling are small, in general agree- 
ment with the values obtained from the early-phase modelling, but indicate the presence of 
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a central low- velocity zone containing ~ 1 M & of mostly carbon and oxygen. The abundance 
of calcium is also very low in the nebula. The pr esence of a dense, ox ygen- rich core has been 
deduced from t he nebular spectra (SN 1998bw, iMazzafi et al.ll200ll ). or via the analysis of 
the light curves jMaeda et alil2003h in all well-observed broad-lined SNe Ic, and it is a strong 
indication that the explosion was not spherically symmetric. Only if the mas s of this zone 



is lar ge can an effect on the early light curve be appreciated (e.g., SN 1997ef, iMaeda et al. 



20031 ). Unfortunately, in the case of SN2006aj spectroscopic observations in the early phase 
extended only to 2006 March 10 (20 days after the explosion), and the evolution of the light 
curve at epochs of 2-3 months, when the linear decline is a signature of the inner dense 
zone (IMaeda et ajj 120031 ). was missed. Nevertheless, the power of a combined approach of 
studying both the early phase and the late phase is highlighted by the present results. 

The nebular lines have rather symmetric profiles. In particular, the [O i] line is not 
as sharp as in SN 1998bw, indicating a low degree of asymmetry. Early radio observa- 
tions suggested that an y aspherical outflow was characterized by a large opening angle 
( iSoderberg et al.1 l2006bl ). and our results co nfirm that the ejecta of SN2006ai look more 
like those of ordinary SNe Ic (e.g., SN 19941, iFilippenko et al.lll995c ISauer et al.l 12009 ) than 
those of SN 1998bw. Early-phase modelling suggested that only a small fraction of the 56 Ni 
synthesised by the SN was ejected at high velocities, and our nebular models confirm that 
most of the 56 Ni was ejected at low velocities. 

The only possible deviation from sphericity may be in the deepest regions, as indi- 
cated by the low- velocity enhancement of the [Ol] line. This may indicate an oxygen- 
rich disk viewed nearly face-on, as in SN 1998bw but of much smaller size, confirming that 



all SNe Ic are prob ably intrinsically aspherical in their cores ([Leonard fc Filippenkol 12005 



Leonard et al.ll2009) . Given the small size of this region, and limited quality of our polariza- 
tion data, we do not expect that we can detect polarization from that structure. Additionally, 
at the late phases of these data, polarization is more sensitive to the distribution of 56 Ni than 
to global asymmetries of the ejecta. The lack of a detection of polarization may therefore 
also indicate that the 56 Ni distribution in the inner ejecta was not very aspherical. 

Our results confirm that the mass ejected in the SN explosion was small, and thus lend 
credence to our earlier conclusion that the progenitor star had a ZAMS mass of ~ 20 M , 
and that probably the remnant was a neutron star (Mazzali et al. 2006b; see also Maeda 
et al. 2007). They may also shed some light on the magnetar phenomenon. If the birth 
of a magnetar was at the origin of the XRF, this magnetic activity did not lead to gross 
asphericity in the SN ejecta or to the synthesis of much 56 Ni at high velocities. Probably, it 
did not lead to a jet at all, but rather to broad polar outflows. This should be useful not 
only for further modelling, but also to understand better the nature of XRFs. As mentioned 
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by iPian et al.l (120061 ). XRFs may be intrinsically weaker and more spherical versions of the 



classical, highly collimated cosmological GRBs. 

Continued monitoring of SN2006aj is planned at both Keck and VLT. When more data 
are collected, we will study in detail the evolution of the spectra and the light curve in the 
late phase. 
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the Explosion of Massive Stars." This work was partially conducted by the groups of J.S.B. 
and A.V.F. under US Department of Energy SciDAC grant DE-FC02-06ER41453. A.V.F. 
is also grateful for the support of NSF grant AST-0607485. 
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Fig. 1.— VLT-F0RS1 spectrum of SN2006aj on 2006 September 19-20. The middle solid 
curve shows the observed spectrum. The upper thin dotted curve is the observed spectrum 
re-calibrated according to the host-galaxy photometry (see Section 2). The host-galaxy 
spectrum is shown as a dashed curve. The spectrum of the SN, obtained by subtracting the 
galaxy spectrum (including the narrow emission lines) from the re-calibrated spectrum and 
binned using a 20 A boxcar, is the lower solid curve. 
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Fig. 2. — The spectrum of SN 2006aj obtained at the Keck Observatory on 2006 July 26, cor- 
responding to a rest-frame epoch of 153 days after XRF 060218 (solid/blue line), calibrated 
and binned with a 25 A boxcar, compared to the spectrum of SN 1998bw spectrum obtained 
on 1998 September 12 (Patat et al. 2001), at an epoch of 139 days after the explosion, 
which is assumed to have occurred on 25 April 1998, in coincidence with the GRB (grey, 
dashed/red line). The flux of SN 1998bw has been rescaled to match that of SN2006aj at 
the peak of the [O i] line. See the electronic paper for a color version of this figure. 
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Fig. 3.— A comparison of the VLT 2006 Sept 19-20 (206 rest-frame days after XRF 060218) 
spectrum of SN2006aj, calibrated and binned with a 20 A boxcar (solid/blue line), with the 
SN 1998bw spectrum obtained on 1998 November 26 (Patat et al. 2001), 214 rest-frame days 
after GRB 980425 (grey, dashed/red line). See the electronic paper for a color version of this 
figure. 
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Fig. 4.— A com parison of the [O i] AA6300, 6363 line of SNe 2006aj (solid/blue line), 1 998bw 
(dotted/red line. lPatat et al.ll200ll ). and 2003jd (dashed/black line. lMazzali et al.ll2005l ). The 
profile of the line in SN 2006aj is less peaked than that of SN 1998bw, although the average 
expansion velocity of other elements is larger in SN 1998bw than in SN 2006aj, indicating a 
smaller degree of asphericity. See the electronic paper for a color version of this figure. 
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Fig. 5. — The Keck spectrum of 2006 July 26 (grey/blue line) and the corresponding synthetic 
spectrum (black line). See the electronic paper for a color version of this figure. 
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Fig. 6. — The VLT spectrum of 2006 September 19-20 (grey/blue line) and the corresponding 
synthetic spectrum (black line). See the electronic paper for a color version of this figure. 
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Table 1. Model Properties 



UT Date 


Epoch 


V 


M( 56 Ni) 




M(O) 


M(C) 


L 




(rest-frame days) 


(kms -1 ) 


(Mq) 


(M Q ) 


(M ) 




(erg s" 1 ) 


2006 July 26 


153 


8000 


0.20 ±0.01 


2.07 ±0.20 


1.50 ±0.15 


0.30 ±0.05 


1.3 x 10 41 


2006 Sept 19-20 


206 


7400 


0.19 ± 0.01 


1.94 ±0.25 


1.42 ±0.15 


0.25 ± 0.10 


5.1 x 10 40 



